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Abstract 

We investigate liquid layers adsorbed at spherical and corrugated cylindrical 
substrates. The effective Hamiltonians for the liquid-gas interfaces fluctuat- 
ing in the presence of such curved substrates are derived via the mean-field 
density functional theory. Their structure is compared with the Helfrich 
Hamiltonian which is parametrized by the bending and Gaussian rigidity co- 
efficients. For long-ranged interparticle interactions of van der Waals type 
these coefficients turn out to be non-universal functions of interfacial curva- 
tures; their form varies from one interface to another. We discuss implications 
of the structure of these functions on the effective Hamiltonian. 

Key words: wetting, Helfrich Hamiltonian, bending and Gaussian rigidity 
coefficients 
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1. Introduction 

The morphology of wetting_ layers adsorbed at planar substrates is nowa- 
days rather well understood @, 0, 0]. However, adsorption at non-planar 
substrates is rather ubiquitous in nature and one is led to study the behavior 
of inhomogeneous fluids in the presence of curved substrates, for example 
the adsorption at large colloidal particles 0, @, 0] and fibers j§, 0, [h| , capil- 
lary condensation in porous media, liquid bridge formations between spheres 



111 . [12[ and cylinders [13j, [14| . Though 60 years passed since the landmark 
work of Tolman [l5| on the surface tension of curved interfaces there are still 
many open problems related to adsorption at curved substrates. 
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Two basic approaches to the problem of the surface tension coefficient 



dependence on interfacial curvature have been developed [16(. The first is 
based on the analysis of the interface fluctuating around the planar config- 
uration. The second approach focuses on the average shapes of the fluid 
interfaces curved around the non-planar substrates. In both cases the sur- 
face free-energy density is expanded in powers of curvature. The zeroth 
order term represents the surface tension coefficient of a planar interface, 
the coefficient in front of the linear term is related to the so-called Tolman 
length 15] while the second order terms contain the bending and Gaussian 
rigidity coefficients. These two coefficients for membranes parametrize the 
phenomenological Helfrich Hamiltonian 17j. One of the still open problems 
in this field is concerned with the status of this expansion and the existence 
of the bending and Gaussian rigidity coefficients. In particular one would like 
to know under what circumstances, if any, one might expect the non-analytic 
dependence of the surface free-energy density on the interfacial curvature. 

In this paper we discuss these issues while investigating the fluctuat- 
ing, cylindrically shaped liquid-gas interface. First (Section [3]) we follow the 
standard approach and calculate the sum of bending and Gaussian rigidity 
coefficients for the spherical geometry. Then we investigate the cylindrical, 
fluctuating interface (Section [5]) and obtain information separately on the 
bending and on the Gaussian rigidity coefficients. We discuss the structure 
of the ensuing coefficients from the point of view of their dependence on the 
interparticle interaction and the substrate geometries ( Section EJ). In the last, 
Section [7J we summarize our results. 



2. The model 

In order to calculate the effective interface Hamiltonian for an interface 
separating a liquid-like layer adsorbed on the curved substrate from the gas 
phase we employ the density functional theory (DFT). The grand canonical 
density functional f2([p(r)], T, //) depends parametrically on temperature T, 
chemical potential fi, interparticle potential w(r) assumed to be spherically 
symmetric, and on the external potential V(r) representing the effect of the 
substrate. The interparticle potential w(r) is split into the short-ranged 
repulsive part Wh s (r) and the long-ranged attractive part w(r) which is of 
the van der Waals type 

w(r) = w hs (r) + w(r) . (1) 
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The mean-field version of the density functional takes the form [18|, [19 
n([p(r)],T,p) = jd 3 rf hs (p(r)) + \ jd 3 rjd 3 r'w(\r - r'\)p(v)p(v') 
+ fd 3 r{V ext {r)-p)p{r), 



(2) 



where the first term on the rhs represents the free energy of the fluid interact- 
ing via short-ranged repulsive potential Whs{r) evaluated in the local density 
approximation. In the following analysis the long-ranged attractive part of 
the potential w(r) is modeled by 

w{r) = 'WT^Y 3 (3) 

which decays ~ r~ 6 at large distances (we neglect the retardation effect). 
The parameter k corresponds to the hard core radius of the fluid particles 
and the amplitude A > characterizes the strength of the attraction. The 
external potential V ext {r) acting on the fluid particle located at position r 
comes from its interactions with all substrate particles. This interaction can 
be again split into the short- and long-ranged part. The long-ranged part is 
again modeled by the potential w s {r) = — A,/(/t 2 + r 2 ) 3 , i.e. 

V ext (r) = [ dr'p s w s (\r-r'\) , (4) 

JVs 

where V s denotes the region occupied by the homogeneous substrate with 
density p s . The effect of the short-ranged repulsive part of the substrate- 
fluid interaction is to prevent the fluid particles from penetrating the region 
V s ; it is taken into account by the appropriate specification of the different 
regions integration present in the density functional, Eq. (J2j). 

The thermodynamic state of the fluid is taken to be slightly off the bulk 
coexistence line in the regime of the stable bulk gas phase, and far away from 
the critical point. This implies that the bulk correlation length is comparable 
with the size of fluid particle characterized by k. In such circumstances the 
position of the liquid-gas interface is represented by function z = /(R), 
R = (x, y), and the nonuniform density profile p(r) can be described by the 
so called sharp-kink approximation 

p shk (R, z) = Pl 9(/(R) - z) + Pg Q(z - /(R)) , (5) 
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where Q(z) is the Heaviside function while pi and p g denote the densities of 
the coexisting bulk liquid and gas phases, respectively. For a finite system 
the density functional fl([/?(r)], T, p), Eq. (jS}, evaluated at p(r) = p a hk(&i z ) 
can be represented as the sum of bulk, surface, line, etc. contributions [19] . 
We discuss two shapes of substrates which induce the corresponding shapes 
of interfaces: a spherical substrate of radius R s , and an axially symmetric 
substrate represented by a corrugated cylinder of infinite length. In each 
case the gas phase is present away from the substrate while a liquid-like 
layer is adsorbed on the substrate. Our analysis is concentrated on surface 
and interfacial contributions to Q([p(r)],T, p). 



3. Spherical interface 

The system under consideration in this section consists of the spherical 
substrate of radius R s on which a uniform liquid layer of constant thickness 
i is adsorbed. The remaining part of the system is filled with the gas phase 
of volume V g while the volume of the liquid layer is denoted by V\. The 
density functional, Eq. (j2j), evaluated for this geometry and for the spherically 
symmetric density profile in the sharp-kink approximation - after taking into 
account contributions due to the finite size of the system - equals (the upper 
index (s) refers to the spherical case and the lower index s refers to the 
surface properties) 

fiW = (V g + Vi)u b (p g )+n^ , (6) 

where u>b(p) is the grand canonical potential density of a fluid system with 
the uniform density p 

Ub(p) = fhs(p) + ~ p 2 Jdrw(\r\) - pp . (7) 
The surface contribution has the following form 

n« = i*#M B) = ^R 2 S + ^(R s ) + ^^^{Rs + £) + 4 s n \(Rs, *)) • 

(8) 

The first term on the rhs in Eq. (jSJ) is proportional to the volume of adsorbed 
liquid layer of thickness £ 

(S) (Rs + () 3 — Rl , x 

ex = m 2 b ^ pl ' ' ( ' 
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The next two terms crj 5 (R s ) and a\ s g (R s + £) are the surface tension coef- 
ficients of the substrate-liquid and liquid-gas spherical interfaces with radii 
R s and R s + £, respectively. The last term describes the interaction between 
the substrate-liquid and liquid-gas interfaces 

w£>(M) = Ap(p^W(a,£,H)-p s ^!(a ) ^H)) , (10) 
where Ap = pi — p g , and 

£2M,M) = ^-J d 3 r / dy«;(|r-r'|) . (11) 

A-Kd z J\ r \> a+e J\ r \< a 

Within the sharp-kink approximation employed in this paper each of the 
two contributions to the interaction between the substrate-liquid and liquid- 
gas interfaces, Eq. fllQI) . factorizes into the term depending on the densities 
only, and the term determined exclusively by the interparticle interaction 
and geometry; the latter one is denoted with a hat. With its help one can 
express all the interface-interface interactions and surface tension coefficients 
in the problem at hand. For example, the expressions for the surface tension 
coefficients take the following form 

o-SV) = -^(p^2(a,0,H)-2p s ^l(a,0,H)) , (12) 



<#V) = -i(Ap) 2 ^l(a,0,M). (13) 



'ig 

In the limit of large sphere radii one obtains the well known expression [ij 
for the surface tension coefficient for the planar liquid-gas interface 



1 f°° 1 

ai;\a - oo) = a® = --(Ap) 2 vr J dr r 3 W (r) = --(Ap) 2 ^ . (14) 

For the van der Waals interactions, Eqs (131T4"|) one gets 0, 



a£>(a, £, H) = ^77^ { 4a(a + £) + 4a 3 arctan ■ 



Qa 2 { v 2a + 

2£(3a 2 + 3a! + I 2 ) arctan ■ 



1 + £(2a + £) 
+ ln(l+£ 2 )-ln[l + (2a + £) 2 ]} , 



(15) 
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where parameters a = a/n, t = 1/ k are dimensionless, and 



(p) 
ig 



Air 
'4^2 



< . 



(16) 



Similarly, the interaction between the spherical substrate-liquid and the 
spherical liquid- gas interfaces tends in the limit 1/a — > 0, I fixed to the 
interaction between two planar interfaces 



[l -£arctan(l/£)] 



(17) 



The approach of uj^ t (a,£,[w]) to Cj^ t {i, [w]) is described by the following 
asymptotic expression 



H) 



1 + - + ... 

a 



;is) 



Using Eq. (1131) for the van der Waals interactions one gets the following 
expression for the liquid-gas interfacial tension coefficient of spherical inter- 
face with radius a 



a 



(p)£ 
ig 3 



1 + a arctan(l/2a) — 



ln(l + 4a 2 
4a~ 2 



(19) 



In the limit of small curvature 1/a the above expression takes the asymptotic 
form 



a 



ip) 
ig 



1 + 4 (-1- 121n2 + 121n~ ] 4 + 



(20) 



The interfacial tension coefficient (a) is a non-analytic function of dimen- 



sionless curvature 1/a. It is smaller than in the planar case and attains 



the planar limit from below. The first curvature dependent correction to the 
planar surface tension coefficient is proportional to ln(l/a)/a 2 . Qualitatively 
similar results have been obtained also for a different model of a long-ranged 
interparticle potential 0, . 



4. Local Hamiltonian 

Fluctuating membranes can be analyzed with the help of a local Hamil- 
tonian proposed by Helfrich 17| . It represents the energy of a membrane 
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in terms of the mean H = (l/i?i + l/R^)/^ and Gaussian G = I/R1R2 
curvatures, where R\ and R2 denote the principal radii of local curvature. 
Hamiltonian of similar form has been also derived for fluctuating liquid-gas 
interfaces 19 



Helf 



J i 2 s [a^ + ku(H — Co) 2 + k G G 



(21) 



where d s denotes the surface element, Cq - the spontaneous curvature, and 
kn, kc are the coefficients of bending and Gaussian rigidity, respectively. 



The parameter Co measures the profile assymetry [2l|, |22j, [23] and for the 
sharp-kink approximation used throughout this paper it is zero. The first 
non-vanishing corrections to the planar interface Hamiltonian are quadratic 
in inverse curvature radii. Accordingly, for the spherical liquid-gas interface 
with radius a the integrand in Eq. (1211) can be written as 



<tJ J + H 2 + k^ G 



"G 



(22) 



where H 2 = G = 1/a 2 . Comparison of the above equation with Eq. (l20p 
shows that the sum of bending and Gaussian rigidity coefficients depends on 
the curvature radius 



4\a) + k 



g y 



(p) 2 l 
19 36 



-1 - 121n2 + 12 In 



(23) 



The logarithmic dependence of the sum in Eq. (1231) on the radius of curvature 
reflects the presence of van der Waals interparticle interactions in the system 
under consideration Id, l20fl . 



5. Fluctuating cylindrical interface 

In this section we analyze the liquid-gas interface which fluctuates around 
its cylindrical configuration. This fluctuating interface separates the liquid 
layer adsorbed at the corrugated cylindrical substrate from the gas phase. For 
simplicity we assume that both the fluctuating liquid-gas interface and the 
corrugated substrate surface have axial symmetry and - in cylindrical coordi- 
nates r_i_, 0, z - their positions are given by r± = f(z), and r± = s(z), respec- 
tively; see Fig.[TJ We aim at deriving the effective Hamiltonian J^rJf] for 
this fluctuating liquid-gas interface. After evaluating the density functional 
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Q([p(r)],T, /i) for the sharp-kink density profile one extracts the surface con- 
tribution 

Qs ([s], [f],T,n) to Q([p(r)],T, fi). It is a functional of both the substrate 
surface shape and the shape of the fluctuating interface (the dependence on 
thermodynamic parameters is not displayed) 




Figure 1: Section of the system with axially symmetric substrate r± — s(z) and the 
liquid-gas interface rj_ = f(z). 



nH*]>\f\) = + <([/]) + ^M,[f\) . (24) 

The first two terms correspond to the free-energy functionals of the corru- 
gated solid-liquid and the liquid-gas interfaces, respectively. The last term 
denotes the interaction between these two interfaces 

^M,[f\) = Ap(pAS(M,[/],H)-pAS(M,[/],H)) , (25) 

where 

AS(M,[fl.M) = fa fa' [ dV / d 2 r>(|r-r'|) . (26) 

J J J\r ± \<s(z) J\r' ± \>f(z>) 

The two-dimensional vectors rj_ and r' ± are perpendicular to the z axis, 
r = (r_i_, z). The range of integration over z and z' is determined by the size 
of the system along the symmetry axis, say [— L, L]. However, we are not 
interested in the finite size effects due to this cut-off along the z axis and 
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extend the integration limits L — > oo keeping in mind, that only the quantity 
^innN' [f\)l L ls we ^ defined in this limit. 

Assuming small corrugation <C 1 one can perform the gradient 

expansion of f2^([s], [/], [w]) up to bilinear terms (Appendix |A|) and obtain 
the local form of the interaction energy 

n&(M, [/], M) = 2vr fdzs(z){u^(s(z),f(z)) 

J (27) 

+ *«(.«,/«)} ■ 

The first term on the rhs describes the interaction between two coaxial, 
undulated cylindrical interfaces in the Derjaguin approximation 24|. The 
function u^ t (a, b) itself describes the interaction of cylindrical interfaces of 
constant radii a and b (a <b), i.e. without any undulation. The second term 
is the correction due to undulation; it vanishes for ideally cylindrical shape 
s(z) = const. This term turns out to be crucial in determining the coefficient 
of surface tension. 

The surface free energy of corrugated, axially symmetric substrate-liquid 
interface with ^-dependent radius r_|_ = s(z) is equal 

)( c )rroA _ _1 ( n . o( c )/Toi m r„„n_o„ o( c ), 



nS'(M) = -2^^^(M,M,M)-2p s ^(M,M,K])J , (28) 

while the corresponding expression for the surface free energy of liquid-gas 
interface of radius r± = f(z) is equal 

<%]) = -^( A P) 2 ^int([f], If], M) • (29) 

For small undulations |/'(-2)| <C 1 one can rewrite [/]>[/]> M) i n the 
local form 

<''° C ([/]) = 2 7 r/d,/(,){a;;)( / (z)) + ^a;; ) (/(z))l , (30) 



where 



^(a) = &&(a, a) = 1 (i?(-4a 2 ) - iT(-4a 2 
a;;)(a) = ^(a,a) = ^2af^. 



(31) 
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The symbols K{x) and E{x) in Eq. (j31j) denote the complete elliptic integrals 
of first and second kind [25j], respectively. The expression a\^{a) represents 
the surface free energy of cylindrical surface of radius a and is often shortly 
denoted as the a-dependent coefficient of surface tension. It is an increasing 
function of radius a, see Fig. [2j On the other hand, the coefficient in front 
of ~/'(z) 2 in Eq. (l30|) is different from <rj£ and is denoted as (a); it is not 
a monotonous function of a, see Fig. [2J It has a single maximum and for 



A 




Figure 2: The radius of curvature dependence of the functions <5yJ (a) (solid line) and 
affi (a) (dotted line) in the a$ units. 



a — > oo it tends - similarly to <jf g (a) - to the coefficient of the surface tension 
of planar interface 



to 



1 + ^ (l-61n2 + 21n-l 4 



(7 



is 



a I er 



1 + 4 ( -3 + 61n2-21ni ) 4 + 
lo \ a I a 1 



(32) 



We recall that the local capillary-wave Hamiltonian for an interface fluc- 
tuating around a planar configuration and whose position is described by 
/ = f(x, y) is - in the limit of small interfacial undulations - given by JJJ, [26 



jdxjdya^ (l + i(V/) 2 ) = jdha^ , (33) 
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where d 2 s denotes the area element of the interface. One might expect 
that the local form of the Hamiltonian for an interface fluctuating around 
a cylindrical configuration will display the structure similar to the one in 
Eq. (1331) . i.e. the coefficients in front of 1 and |(V/) 2 in Eq. (13U1) will be the 



<r\ftf) in Eq. ([33]). However, 



same which would amount to replacing 

this is not the case because af g \f) ^ (/), and the Eq. (1311 cannot be 
rewritten in a form similar to Eq. (1331) . 

To discuss this issue from a somewhat different perspective we return 
to the Helfrich Hamiltonian (Eq. (1211) ). For an axially symmetric interface 
rj_ = f(z), the mean and Gaussian curvatures are given by 




1 



G 



r 

/(i + n 2 ' 



(34) 



and the Helfrich Hamiltonian expanded up to terms quadratic in f'(z) - after 
some algebra, integration by parts - is equal to 



u(<0 
"H 



+ I ~ 



4/2 1 \"ig 4/2 j 2 

+ 7TAir +k ° 



(35) 



(c) (c) 

where the last term indicates that we allow the rigidity coefficients k H , k G 
to be considered as functions of the local interface position f(z). Comparing 
the above expression with the local form of the surface free energy obtained 
within the density functional approach, Eq. ( |30|) . one gets the following ex- 
pressions for the bending and Gaussian rigidity coefficients 



k 



h y 



k£(a) = -a? g >K'-\n- + k 



ct^k 2 ^ (l -6hi2 + 21n- 



(36) 



Note that the J^^ lf [f] depend ence on the Gaussian rigidity coefficient, 
Eq. (1351) . is only via its derivative and thus the function k G is known only 
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up to an arbitrary integration constant k$. The value of this constant is ir- 
relevant because the integral of the Gaussian curvature over entire interfacial 
area is zero on the basis of the Gauss-Bonnet theorem j2?| . The logarithmic 
terms in rigidity coefficients arise from the corresponding logarithmic terms 
in surface tension coefficients in Eq. fl32|) . 



6. Bending and Gaussian rigidity coefficients 

For short-ranged forces the mean and Gaussian rigidity coefficients which 
appear in the Helfrich Hamiltonian (Eq. (I2TI) ) are constant. They do not 
depend on the actual interfacial geometry, and are usually derived either by 
considering the thermal fluctuations around the planar configuration, or by 
investigating the mean shapes of interfaces curved around cylindrical and 
spherical substrates. In the first approach based on the density functional 
theory supplemented by the sharp-kink approximation these coefficients are 
obtained as the fifth moments of the attractive part of the interparticle in- 



teraction 19 



POO 

h$ = — (Ap) 2 y drr 5 w(r) 



4 P) = (37) 

while the surface tension coefficient is represented by the third moment of 
this interaction 



°19 



(P) 



IT f°° 

--(Ap) 2 / drr 3 w(r) . (38) 
2 Jo 



The second approach consists of two steps. First, from analyzing the cylin- 
drical shape of the liquid-gas interface (with zero Gaussian curvature) one 
derives the expression for the bending rigidity coefficient. Second, from con- 
sidering the spherical interface (Section|3j) one gets the expression for the 
sum of bending and Gaussian rigidity coefficients. Assuming the geometry 
independence of these coefficients, one can obtain also the expression for the 



Gaussian rigidity coefficient [28|, |29|, [30, [3JJ . 

The bending and Gaussian rigidity coefficients derived in this way are 
equal to those derived in the former approach, provided the fifth moment 
of attractive part of the potential exists. On the other hand in the case of 
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van der Waals interactions the rigidity coefficients for cylindrical and spher- 
ical geometries turn into functions of the mean curvature of the interface 
This fact has been already noted in the literature. 



29. 30 



However, 

the bending and Gaussian rigidity coefficient functions were derived within 
the second approach 29, 3^] by assuming that they depend on the interfacial 
curvature in a universal way, i.e. in the same way for the cylindrical and 
spherical geometry. 

By investigating the liquid-gas interface fluctuating around the corru- 
gated cylindrical substrate we derived expressions for both the bending and 
Gaussian rigidity coefficients. The sum of these quantities differs from those 
obtained by considering the spherical interface. It means that for van der 
Waals forces the bending and Gaussian rigidity coefficients are geometry 
dependent. We stress that these conclusions apply only for interparticle po- 
tentials for which the cylindrical Gaussian rigidity coefficient depends on the 
curvature of the interface, e.g. for potential in Eq. (j5J. Otherwise the Gaus- 
sian rigidity coefficient cannot be derived using this method because of the 
Gauss-Bonnet theorem. 

The above non-universal behavior of bending and Gaussian rigidity coef- 
ficients is restricted to the type of interparticle potential considered in this 
paper. Analogous considerations on the level of the liquid-gas surface tension 
coefficient for interparticle potentials with r -4 decay at large distances, e.g. 
w(r) = —A/ (K 2 + r 2 ) 2 , lead to the conclusion that already the surface tension 
coefficient in the planar geometry (Eq. (1551) ) does not exist. The surface ten- 
sion coefficients in cylindrical (cr^(i2 c )) and spherical (a\ S g(R s )) geometries 
can be expressed in terms of elliptic integrals and tend to infinity in the limit 
R C)S — > oo. The dominant terms take the following form 



— (Ap) 2 h -31og2 + log — 

Air , A . n ( , k 

- — (Ap) 2 f-log2 + log — 



R c — »oo 



R s — >oo 







(39) 



Thus also in this case one cannot define the curvature dependent surface 
tension coefficient which is universal in this sense that it does not distinguish 
from which geometry the curvature comes from. 
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7. Summary 

The effective Hamiltonian for the liquid-gas interface fluctuating around 
the corrugated cylindrical substrate has been derived within the mean-field 
version of density functional theory for fluid whose particles interact via van 
der Waals potential. This Hamiltonian is non-local but in the limit of small 
interfacial undulations and small cylindrical substrate corrugation it takes 
local form. However, the structure of the effective Hamiltonian, Eq. fl30|) . in 
its local form is qualitatively different from the one describing the fluctuating 
interface around the planar configuration Eq. fl5H|) ; it cannot be presented as 
the integral of the local curvature dependent surface tension coefficient over 
the interfacial area. 

The local effective Hamiltonian for the fluctuating cylindrical interface 
can be rewritten in the form proposed by Helfrich only when the rigidity 
coefficients multiplying the mean curvature squared and the Gaussian curva- 
ture are replaced by functions of the local curvature. The structure of these 
functions turns out to be non-universal, i.e., their form depends on whether 
one considers spherical or cylindrical interface. 
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A. The gradient expansion of the interaction energy between two 
corrugated cylindrical interfaces 

The interaction energy of two axially symmetric interfaces r± = s(z) and 
r± = f(z), f(z) > s(z) is equal 



«&([*]>[/]>])= Idz Idz' 



d 2 rj_ 



r±\<s(z) 



d 2 r>(|r-r'|) 



'1 !>/(*') 



Jdz Jdz' [jV(s(z) : oo : z- z') - J^(s(z)J(z') : z - z')] , 

(40) 



where 



j( c \a,b,z') = [ d 2 r ± [ d 2 r' ± w ( J(r ± - r' ± ) 2 + z' 2 j 

J0<\r ± \<a J0<\r' ± \<b V / 



(41) 



denotes the interaction between two coaxial discs with radii a and b separated 
by distance z'. Assuming small corrugations, i.e., \f'(z)\, \s'(z)\ <C 1 one 
obtains 

nS[«,/l= J dz Jdz'[J(s(z),oc,z-z')-J(s(z'),f(z),z-z')\ 
= Jdz Jdz' [J(s(z),oo,z')-J(s(z' + z)J(z),z')} 



= dz dz' 



.12 



J(s(z), oo, z) - J{s(z) + s'(z)z + s"(z) — , f(z), z') 



dz / dz' [J(s(z), oo, z') - J(s(z), f(z), z')\ 



s'(z) 2 



d 2 J(s(z)J(z),z') 1 d dJ(s(z),f(z),z'Y 



ds(z) 



s'(z) dz 



ds(z) 



-Jdz Jdz' 

=2«fdz8(z)&&(s(z),M) + fdzfd* 2 - d8 (z)d(f(z)) 
=2vr Jdzs(z) (*(*), /(*))+ AW (*(*),/(*))} • 



/S '(£)f(£) /2 9 2 J( S (^),/(^),z') 
— — z 



(42) 
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In the last line, the first term in the bracket 6j^ t (a, b) denotes the interaction 
between two cylindrical interfaces of radii a and b, respectively 2p| 



An 
4k 2 



dz' 



m 2 + n(z' 2 + 1) 



2a(z' 2 + l) 2 \ ^(z'l + 1)2 + 2n (z' 2 + 1) + 



m 



(43) 



where m = b 2 — a 2 , and n = b 2 + a 2 (b = bj k, a = cl/k). The second term 
X^(a,b) describes the coupling between the corrugations 



A (c) (a,&) 



An 



2bE 



(b-a) 2 +l 



(44) 



In the limit a, b ^> 1, 
values 



4«V ((6 + a) 2 + l)v/(6-a) 2 + l 
6 — a fixed, both quantities tend to its planar 



(45) 
(46) 



where is given by Eq. (117)1 . and 



A (p) (£) 



A7T 



4k 2 y i + £ 2 

On the other hand, for a fixed, and I 3> 1 one has 

A( c )(a,6) = \V>(£) 
and from numerical computation one obtains 



(47) 









1 - 2ai + 


— a - 











^7 + 



(48) 



(49) 
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